Entanglement and area laws in weakly correlated gaussian states 
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We examine the evaluation of entanglement measures in weakly correlated gaussian states. It 
is shown that they can be expressed in terms of the singular values of a particular block of the 
generalized contraction matrix. This result enables to obtain in a simple way asymptotic expressions 
and related area laws for the entanglement entropy of bipartitions in pure states, as well as for 
the logarithmic negativity associated with bipartitions and also pairs of arbitrary subsystems. As 
illustration, we consider different types of contiguous and noncontiguous blocks in two dimensional 
lattices. Exact asymptotic expressions are provided for both first neighbor and full range couplings, 
which lead in the first case to area laws depending on the orientation and separation of the blocks. 

PACS numbers: 03.67.Mn, 03.65.Ud, 05.30.Jp 



I. INTRODUCTION 

Entanglement is a valuable resource that plays a key 
role in quantum information processing and transmission 
based on qubits [H-[l] or on continuous variable systems 
It has also provided new insights into the role of 
quantum correlations in the critical behavior of many- 
body quantum systems [6l-fTojj . Nonetheless, the evalu- 
ation of genuine quantum correlations for a given state 
of a many-body system is in general a difficult task. On 
the one hand, rigorous computable entanglement mea- 
sures exist just for pure states, where the entanglement 
entropy, i.e. the entropy of the reduced state of a sub- 
system, pro vides the basic measure of bipartite entan- 
glement In the case of mixed states, rigorous mea- 
sures like the entanglement of formation [12( , which is the 



convex-roof extension of the previous measure 
volve a minimization over a very high dimensional space 
of parameters and are therefore not directly computable. 
This has turned the attention to the ne gati vity fl4j . or 
cquivalcntly, the logarithmic negativity |14l . Il5| . which 
quantifies the violation of the positive partial transpose 
separability criterion by entangled states and is a bipar- 
tite entanglement monotone 1A\ , computable in principle 
for any bipartition in any pure or mixed state. Neverthe- 
less, the accurate evaluation of these quantities demands 
a deep knowledge of the many body state, which requires 
in general an amount of information which increases ex- 
ponentially with the system size. This fact limits the 
possibility of closed evaluations to states characterized 
by a manageable number of parameters . 

A prime example of such states are the gaussian states, 
i.e., ground or thermal states of stable gapped Hamiltoni- 
ans quadratic in boson operators, or equivalcntly general- 
ized coordinates and momenta (EflflGj}- For such states, 
of crucial importance for continuous variable quantum in- 
formation |5| , the entanglement entropy of bipartitions of 
pure states and the negativity between arbitrary subsys- 
tems in pure or mixed states can be evaluated in terms of 
the elements of the covariance matrix fl8l - [22l ]. i.e., of the 
generalized contraction matrix of pairs of boson operators 



[23l l24j . However, even in this scenario, the extraction 
of analytic expressions for these quantities for arbit rary 
subsystems is in general not straightforward [IH, EH, |25| . 

The aim of this work is to discuss the evaluation of the 
previous measures in weakly correlated gaussian states, 
such as typical ground states of gapped hamiltonians, 
which can be characterized by excitations over a product 
state. Gaussian states are usually described in terms of 
the phase space formalism p|, which allows to connect 
their entanglement properties with correlations in phase 
space. Here we will consider a different approach, based 
on the Fock representation, which provides an equivalent 
yet in many cases more clear way to evaluate and repre- 
sent entanglement measures (24[. We will show that the 
entanglement entropy and negativity can be expressed in 
terms of the singular values of sub-blocks of basic con- 
traction matrices, which can be evaluated analytically in 
the perturbative limit for some typical couplings. The 
formalism also allows the straightforward derivation of 
area laws [T3, EH, [26l - [28j for these quantities. The emer- 
gent area laws for the entanglement entropy and negativ- 
ity are different, i.e., they depend on distinct measures 
of the boundary size, and are affected by the orientation 
and separation of the subsystems. Let us also remark 
that the ground state of weakly interacting spin systems 
can also be described by gaussian states throu gh d ifferent 
approximate bosonization techniques (23l. |2"H. |29| . entail- 
ing that the scope of the present scheme is quite general. 

The formalism is described in section II, while section 
III considers its application to specific systems, essen- 
tially ground states of two dimensional lattices with short 
range couplings, although the full range case is also con- 
sidered. The present scheme allows to easily obtain ex- 
act analytic asymptotic expressions for the entanglement 
entropy and logarithmic negativity of different types of 
bipartitions and block pairs, both contiguous and sepa- 
rated, which will be compared with exact numerical re- 
sults. They clearly show the emergence of precise area 
laws. Conclusions are finally drawn in IV. We also in- 
clude appendices containing the details of the perturba- 
tive expansion for the symplectic eigenvalue problem and 
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the evaluation of singular values. 



II. FORMALISM 



The blocks of the contraction matrix acquire then a 
simple form in terms of U, V and the diagonal block F' + : 



F~ 



VU l - 
VV^ 



Entanglement entropy and negativity in 
Gaussian states 



The class of Gaussian states in a bosonic system can 
be defined as those states of the form 

P = Tr C , P \- fm T (*) exp(-/3ff)Tt(a) , (1) 

where H is a positive definite quadratic form on boson 
operators b h b\ ([6», &t] = S lj7 [b^bj] = 0), 

H = J2(\i6ij - A+X&J&j + - ±(At.&^ + Ayb}b\) 
A - A+ -A- 



For a pure state, F' 
F+ 



VV\ implying 

F~F~ 



VF^U 1 + UF' + V f 
VF' + V^ + UF' + U^ . 
= and Eqs. <^ lead to F~ 



F 



(6a) 
(6b) 
VU\ 

(7) 



For such states, the entanglement between any subsys- 
tem A and its complement A can be measured through 
the von Newmann entropy of any of the reduced states: 



£ 



A,A 



S{pa) = S(px), 



(8) 



-A- A-A+ 



(2) 



where S(p) = — Trplogp. In a Gaussian state, the va- 
lidity of Wick's theorem[23 implies that the state of any 
subsystem A is also Gaussian and hence fully character- 
ized by the corresponding contraction matrix T>a, which 
is just the sub-block of V with indexes belonging to A: 



with Z = (£ t ), and T(a) = fli ex P(<5A - aj&J) is a 
displacement operator {T (a)biT^ (a) = b{ + oti). In @, 
A is the diagonal matrix of local bare energies and 
Aj • are the coupling strengths between pairs of different 



The von Newmann entropy © can then be expressed in 
terms of the symplcctic eigenvalues f£ of T>_\ as 



F 



A+, At 



AT). In the pure state limit S{p A ) = ^ h(f£ ) , h(x) = -xlogx+(l+x)log(l+x) , 



bosons (A^ 

f3 -> oo, p -> T(a)|0)(0|Tt(a), with |0) the ground state 
of i7. The displacements a, can be taken into account 
by local shifts b{ — > b t — oti, so that in what follows we 
will set on = 0, such that (bi) p = Tr pbi = V i. 

The key property of these states is that by means 
of Wick's theorem (23|, the expectation value of any 
bosonic operator (and hence p) is fully determined by 
the displacements and the generalized contraction 
matrix 



(10) 

In the case of a mixed state or for pairs of non- 
complementary subsystems B, C, the subsystem entropy 
is no longer a measure of quantum correlations. Instead, 
a well known computable entanglement monotone for 
such systems is the negativity N^fi [01 j which is just the 
sum of the negative eigenvalues of the partial transpose 
PbC' associated quantity is the logarithmic negativity 



= iog(i + m 



V=(ZZ^) -M = 



F? 



(b}b 



FT 



F^ 

F~ 
(bibj 



F~_ 
1 + F+ 



B,C) 



log||Pe B c 



(11) 



(3) 



where M = ZZ^ - [(Z^Z*]* = (J_5) is the symplectic 



FZ, F: 



F 



metric and F± -, 4 - 

We may diagonalizc T> or n by means of a symplcc- 
tic transformation W satisfying W^AiW = M, corre- 
sponding to a Bogoliubov transformation Z = WZ' to 
boson operators Z' = (*' t ), such that V = WV'W\ 

with V diagonal (F'+ a , = f a 5 aa >, F'~ = 0). This leads 
to the standard diagonalization of the matrix DAi (as 
\VVMW~ 1 = V'M). The matrix W can be written in 
block form as 



where \\A\\i = tr V A^A denotes the trace norm. For a 
gaussian state, ||pg^||i can be expressed in terms of the 
negative symplectic eigenvalues /„ ,c of the contraction 
matrix T>bc determined by p^ c , with blocks 



r BC — 



F ± F T 
F T 



B C 

F^ 



(12) 



where Fq c denotes the matrix of elements F^ with ieS 



?± 

B 



and j € C, and F t 

E 9{fS fi 



F ± 

B,B ' 



C B.C 



Eq. (jTTJ) then becomes 

-log(l + 2a;). (13) 



9\x 



w = 



U V 
V u 



where U and V should satisfy 



uv t - 



(4) 



(5a) 
(5b) 



We notice that j% c > -1/2 @]. In the case where 
B = A and C = A, Eq. (fT3| can be expressed in terms of 
the symplectic eigenvalues f£ of T>a as 



= 2 E lo §( v-^ + v 1 + /« ) ' ^ 14 ) 

a 

which, like Eq. (J5J), is again a concave function of the f£. 
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B. Weakly correlated pure Gaussian states 

We consider now the case of weakly correlated pure 
Gaussian states, i.e., states in which the local symplectic 
eigenvalues (those corresponding to a single mode A = i) 



(15) 



satisfy fi <C 1 V i, such that each mode is weakly en- 
tangled with the rest of the system. In the local basis 
where = (this implies replacing 6j — > mh — e^Vib], 

with Ui , Vi = \J F " +1 2/f-fi' + 1 ana - the phase of 
F^), fi = F^ and weak coupling implies, together with 
the positivity of V and F+ , that \F^\ < ^/f~f~ < 1, 
< Mm[fi,fj] + fifj < 1 V ij. In this limit, Eqs. 
([7]) then lead, neglecting terms oc (F~F~) 2 , to 



F + w F~F~ 



which for a subsystem A implies 



F A F A 



F~ -F- 
A,A A,A 



(16) 



(17) 



we may neglect its second order effect in Gg and Gc in 
(EH), such that 



F~ F~ 

B,BC BC& 



G, 



F~ F— 

C,BC BCfi ' 



(23) 



depend just on the correlation with the environment of 
BC. If the sites of B and C correlated with each other 
have correlations of the same order (or less) with BC, (i.e. 

II-Fggclloo and H-^cW 00 ° f thC Sam ° ° rdCr aS H^clU' 

at least for the subsets of B and C mutually correlated) 
we can directly neglect Gb and Gc in Eq. ([21]) at order 
H^eclloo- The negative symplectic eigenvalues of T>b,C 
will then be given by minus the singular values of 
Fq C (see Appendix IB"]) : 



■I a 



B,C 



r B,C 



(24) 



which depend again just on the F~ contractions between 
B and C. For instance, this is the case of contiguous 
blocks in a scenario of short range couplings, and also 
that where C is the complement of B (C = B). 

In the general case, however, the whole matrix (|2ip 
should be diagonalized. First order corrections lead to 



Using Eq. (jTT)) and the results of Appendix [3] the sym- 
plectic eigenvalues of T>a will then agree at this order 
with the standard eigenvalues of the matrix 



FaFa 



F~ -F- 
A,A A,A 



(18) 



which are just the square of the singular values cr£ of 
FT t (see Appendix [Bj . We then obtain, at this order, 



A,A\2 



(19) 



Entanglement depends at this level just on the F con- 
tractions between A and A. For instance, in the case of 
a single site i, Eq. ([19]) implies fi m cr 2 , = Yl&i I ^7 1 2 - 

In this regime we may just set h(x) f=s — x\og 2 (x / e) in 
Eq. (p~0|) . such that the entanglement entropy becomes 



--A.A 



-E(^) 2 log[(^)Ve 



(20) 



Considering now the negativity, in the present regime 
the symplectic eigenvalues of T>bc will be given at leading 
order by the eigenvalues of (see Appendix |A"|) 



F BC 



FbqFbc 



Gb F, 



F 



where, for S 



B orC, 
Gs = 



F + 



C.B 



FsFs 



B,C 

G c 



(21) 



(22) 



For pure global states, Eq. (jTHJ) leads to Gs ~ F g ^Fg s , 
indicating that Qs takes into account the correlations 
with the environment of S. Up to first order in Fg C , 



./ a 



B.C 



r B,C 



[(Ge)cra + (Ge)aa]/2 ; 



(25) 



where (G B ) QQ = UiG B U a , {G c ) aa = V^G B V a are the 
diagonal elements in the local basis of B and C where 
(Jg C )oa' = o-®' c S aa > (see App. IB"]). As Gb and Gc are 
positive matrices in the approximations (|22[) ()23[) . neg- 
ative eigenvalues can only arise if (Gb) qq and {Gc) aa 
arc not much larger than erf' C . A sufficient condition 
ensuring a negative eigenvalue f^ c of ([2"Tj) is 



B,C 



> v (Gb)q Q (Gc) 



(26) 



In the present regime, the logarithmic negativity can 
be obtained setting g(x) ~ — 21og(e)x in (]T3"]) . such that 



£# c «-21og(e) J2 f«' C > 



(27) 



i.e., £qc ^ \ \Fb ell 1 m ^ ne approximation (j2"4"]) . For com- 
plementary subsystems (C = B = A), it is verified that 
identity between Eqs. (fbfj) and ([2~T|) holds at leading order 
in the approximation ([19"]) (log(a + \J\ + er 2 ) sa log(e)tr). 



C. Ground state correlation matrix in the weakly 
interacting case 

A particular case of the previous results arises when 
we deal with the ground state of a Hamiltonian of the 
form @. For weak couplings A* <C A, the diagonalizing 
symplectic transformation W such that W^HW = f2©r2, 
with Q aa > = 6 aa 'UJ a , can be evaluated pcrturbatively. At 
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leading order (see Appendix [A} the block U in ((4]) is a 
unitary matrix that diagonalizes A — A + , while 



Vi 



i/3 



(C/tA-[/) Q(3 



(28) 



Note that, in contrast with the conventional perturbation 
theory, a possible degeneracy in the local energies A& will 
not spoil this result if the system is stable (w Q > V a). 
Notice, however, that U can depart considerably from 
the identity if the Xk are degenerate. 

If all local bare energies are nearly equal (|Afc — Xj\ <C 
Afc + Xj « 2A) , and if energy corrections arising from A + 
are neglected (ui a ~ A), Eq. (|28|) reduces to V ~ ^A~U. 
In such a case, Eq. © leads to 



F~ 



~2X 



(29) 



with F + given by Eq. p^|) . In this regime, correlations 
are hence proportional to the pairing couplings A - , de- 
creasing as A^ 1 for increasing local energies. Noteworthy, 
the strength of the hopping interaction A + does not af- 
fect the ground state correlations at this order. When 
non degenerate, it just affects F ± dressing the bare pair- 
ing interactions. 

In the same way, for a common local bare energy A, 
inclusion of second order terms in the couplings leads to 



A" 
~2A~ 



A+A- 



A-A+ 



4A 2 



(30) 



This expression is useful in the present scheme when the 
first order term vanishes (modes i,j unconnected by A - ). 
As the counter terms Gg and Cq in (|2"5l) will be of sec- 
ond order in F~ (Eq. ([23]) ). subsystems unconnected by 
A - but connected at second order through Eq. (|3U|) may 
exhibit an 0(A/A) 2 non-zero negativity if Eq. (|26|) holds. 



\\M AA and Eqs. (fl9" ]) -([2"0 |) lead then to 



£ 



A,A 



\dA\s 



(32) 



at leading order in A /A, which coincides exactly with 
the result in [l9j for non critical harmonic systems. 

The logarithmic negativity presents, however, a 
slightly different behavior: for contiguous subsystems, 
the same procedure and Eqs. (|2~4")) - (l2"7| lead to 



where the boundary measure is now 



(33) 



\dA\i 



Tr 



\M 



A,A\\l 



E7±A,A 



(34) 

the singular values of the matrix M A A (in 
comparison, \\M A ^||| = Ea^o" 4 ) 2 )- ^ eacn s ^ e m A 



with 5 A,A 



has at most one neighbor in A (the opposite may not 
hold) , the rows of M A A will be orthogonal and the sin- 

A , leading to 



gular values will be df'"^ = 



\dA\ 1 = J2J(M. 



l A,A M A ,A 



On — E 



(35) 



i£A 



ieA 



which will differ from (|3"Tj) if n A > 1. In general, Eq. ([35 



may provide a rough approximation to the area Q34[) . In- 
terestingly, in an isotropic hypercubic lattice in d dimen- 
sions, the approximation (|35|) is just proportional to the 
cuclidean area for large planar surfaces, both parallel and 
tilted (with an angle of ir/4 with respect to the principal 
axes of the lattice), which is not true in the tilted case 
neither for |cL4|2 nor \dA\i (see next section). 

In general, for two contiguous subsystems B, C, previ- 
ous expressions generalize to 



D. Area laws 



£^ c cx \dBC)dC\ 



(36) 



The formulation of the area law for systems with local 
interactions starts with the definition of a suitable mea- 
sure for the size of the boundary dA of the subsystem A 
fiol Qjl, HU . An example of such measure is given by the 
number of pairs of first neighbor modes, with one mode 
belonging to A and the other to A. If we define the ma- 
trix M with entries Mij = I if modes i and j are first 
neighbors and otherwise, that measure can be written 
as 



\dA\ 2 = Tn A 



ieA 



= Tr [M AA M AiA ] 



\M 



A,A 



(31) 



where 



= (M A A M A A )u is the number of first neigh- 
bors of mode i in A. For the ground state of a gapped 
bosonic system with constant and isotropic first neigh- 



bor interactions A„ 



4±A% Eq. 



(J2SJ implies F 



AA 



at leading order in A, where 



\8BndCU = \\M BtC \\i=J2^ 



(37) 



is a measure of the contacting area between B and C. 
Again, if each mode in C is linked with at most one mode 
in B, af' C = y^nf, where nf = (M b ,cMc,b)u is the 
number of first neighbors of i in C. 

Previous geometric-like expressions can of course be 
also applied to a general constant coupling A~ = 
iA - Mjj, where Mij = 1 if pairs i,j are linked by 
the coupling and otherwise, leading to effective areas 



|M_4^||i and \dA\ 2 = ||M 



A,A\\2- 



On the other 



hand, they cannot be directly applied to higher order ef- 
fects, like those depending on Eq. (|3U|) . as discussed in 
the next section. 
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FIG. 1. (Color online) The complementary partitions con- 
sidered in Eqs. (|42[) . a) Single site, b) square block parallel 
to the principal lattice axes, c) tilted square block and d) 
checkerboard. 



in the case of a single site (Fig. la) we obtain 



(38) 



In the rectangular n x x n y block lb (parallel to the prin- 
cipal axes), there are three different singular values, cor- 
responding to the horizontal and vertical sides and the 
four corners, given below with their multiplicity: 



T AA\2 



2(n x - 2) 
2(n y - 2) 

4 



(39) 



In the n x n square block lc tilted 45° with respect to the 
principal axes, we obtain, by means of a discrete Fourier 
transform and neglecting corner effects (see Eq. (|C3j) ). 



A.A 



2a x a y cos ^ , 



(40) 



III. EXAMPLES AND ASYMPTOTIC 
EXPRESSIONS 

We will now use the present formalism to obtain an- 
alytic asymptotic expressions for £4^4 and E^^ for typ- 
ical subsystems A, B and C of a two-dimensional lat- 
tice, which will be compared with the exact numerical 
results and the estimations (|32l) - (|36)l . We first consider 
the ground state of a bosonic square lattice with attrac- 
tive first neighbor couplings 



+ 5. 



fj,=x,y 



where tt ;i denotes the unit vector along the p, axis. We 
have considered in Figs. (JTH5|) the isotropic case A^ = 
Ay = A ± , with A _ /A + = 2/3, and a uniform single 
mode energy Xi = A. Away from the critical point A = A c 
(where, for fixed A*, the lowest energy uj a vanishes), the 
system is gapped and a finite correlation length £ < 00 
is expected. Approximately, A c ps 2(A+ + |A~|) (exact 
result for the cyclic case |24|). 



A. Entanglement between complementary 
subsystems 

We first consider the four global (A, A) bipartitions 
depicted in Fig. [U Eqs. (HU), and $Z§§ lead to ana- 
lytic asymptotic expressions for the corresponding singu- 
lar values &a'^- At lowest order, their number is just the 
number of sites at the border. Defining the basic single 
link singular values 

o> = l A u l/(4A) , n = x,y, 



where k = 1, . . . , m and m = 4n — 4 is the number of 
sites at the border. Corner effects will affect essentially 
just 4 of these eigenvalues with O(l) corrections. 

Finally, in the checkerboard partition Id, an exact ana- 
lytic expression for the n x n y /2 singular values is available 
in the cyclic case again by means of a discrete Fourier 
transform (see Eq. (41) in (29j): 



A.A 



J2 a » 



cos ■ 



27rk u 



fi=x,y 



(41) 



,n y /2. 



where k — (k x , k y ) with k x = 1, . . . , n x , k y = 1 

These expressions, together with Eqs. ([19]), ([20]) . (|2"3|) 
and (|2"T|) . allow to easily obtain the asymptotic values of 
the entanglement entropy and negativity of the present 
bipartitions for large A and n. For instance, in the 
isotropic case Aj = A* considered in the figures, set- 
ting a = <7 M and neglecting corner and border effects 
(which just add terms of relative order n _1 ), we obtain 



(42a) 
(42b) 

-(4n)2a 2 log 2 (f2|i), (42c) 



•• 4 e 



for the entanglement entropy of the single mode, the 
parallel and tilted n x n square blocks and the n X n 
checkerboard of Fig. [T] We have replaced sums over k in 

c-d by integrals (££=i/(^) « ff/t^^ Note 
that in the chekcrboard case the entanglement entropy 
scales with the "volume" n 2 of A rather than the "area" 
n, since all links are broken by the partition (maximally 
entangled bipartition (29j|). 

The corresponding values of the scaled logarithmic 



1,1 




C A,A 


i -Ana 2 log 2 ^ , 


C A,A 


■i —8na 2 log 2 a 1 = 


£ d - « 
°A,A 


i-2n 2 \og 2 {a 2 e) = 




0.00 



klX r 



FIG. 2. (Color online) Exact and asymptotic (Eqs. (|%2" ]l -(|l3" )l ) 
results for the scaled entanglement entropy £_4,_4 (top) and 
logarithmic negativity ^ (center and bottom) of the four 
bipartitions a, b, c, d, of Fig. [TJ as functions of the ratio A/A c . 
In the top and central panels results were scaled with the 
boundary measure \dA\2 (A, An, 8n, 2n 2 in a,b,c,d, according 
to Eqs. (|31|) - (|46|l 'l. which is seen to provide an adequate scal- 
ing for £ A ^ but not £ A \. The latter scales accurately with 

the measure |cM|i (2, An, ^n, according to Eqs. (f34|) - 
((46J ) , as verified in the bottom panel. Results correspond 
to a 30 x 30 lattice with A"/A + = 2/3 and 10 x 10 blocks in 
b) and c) . 



negativity £* A = £^/[2 log(e)] are 



£ 



Wb_ 
A,A 



Ana . 
16 



£ 



' AT d 
A,A 



7T 



-~^na=(AnW2^a, 
■2 



- = (y)V4(^) 2 ^. 



(43a) 
(43b) 

(43c) 
(43d) 



The last expressions in (|42 |) -(|43 |) indicate the way to 
read them. They are of the form 



£ A ,A~-Lma 2 log(a^) 



(44) 
(45) 



where L is the number of modes at the border of A (L = 
1, An, An, n 2 /2), m is the number of connections with the 
environment A per mode at the border (in = 4, 1,2,4), 
i.e., the number of links per mode broken by the parti- 
tion, and a j , f3 j , with a = e/2 « 1.36, /3 = 2^/2/it w 0.9, 



FIG. 3. (Color online) Non complementary subsystems: Con- 
tiguous (top) and one-site separated (bottom) blocks, with 
contacting sides parallel (left) and tilted (right) with respect 
to the principal axes. The negativity and its size dependence 
are determined by both separation and slope of the contacting 
boundary. 



are geometric correction factors for the tilted (j = 1) and 
checkerboard (j = 2) cases (j = for the single mode and 
parallel square). We can easily identify from (|44p - (|4"5")) 
the boundary measures of Eqs. (|32l) - (|33l) : 



\dA\ 2 = Lm, \dA\i =L^ip j 



(46) 



Comparison with the exact numerical results (Fig. [2]) 
indicate that all these asymptotic expressions are actu- 
ally quite accurate already for A > 4A C . The scaling of 
£ A x with the area \dA\i rather than \dA\2 is clearly ver- 
ified. Moreover, this scaling is more accurate than that 
of the entanglement entropy £ A ^ with |cL4.|2, since the 
latter contains in lc — Id an additional geometric correc- 
tion Lm\o%{a? 'm)a 2 (Eq. (|4"4")) ). not comprised in (|32p . 
Note also that in the case of the tilted block, \dA\2 = 2L 
and \dA\i = L\/2(3 = (A/-k)L are, respectively, larger 
and smaller (90%) than the geometric perimeter \J~2L. 

We may also rapidly determine with Eqs. (|3"9")l and 
(|4"4" f -- (|4"5"|) the corner effects in case lb. The actual asymp- 
totic expressions for the finite n x n parallel block are 



£ b - 

C A,A 



£ 



■Afb_ 
A,A 



-4(n-l)a 2 log 2 ^-4^1og 2 ^, (47) 



(48) 



A(n - l)a + 4(V2 - l)cr 



where the first term is proportional to the number of sites 
at the border, A(n— 1), and the second represents the pos- 
itive correction arising from the four corners, reflecting 
their increased coupling with the environment A. 



B. Non-complementary subsystems 

Let us now consider the non-complcmcntary subsys- 
tems of Fig. |31 For contiguous parallel blocks contacting 
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at n sites, F B c has n identical singular values 

o** = <r*. (49) 

In the case of contiguous blocks with contacting surfaces 
tilted 45° with respect to the principal axes, we obtain, 
neglecting edge effects, the same expression ([4*U|) for the 
cr, ' , with m replaced by the number of contacting sites 
n and k = 1, . . . , n. In the isotropic case we then obtain, 
using Eqs. flMHUD; 

gggfantr, (50) 

for the logarithmic negativity of parallel and tilted con- 
tiguous blocks, which are clearly of the form (|45j) or (|33l) : 
|9Kn9C|i = n and 4n/7r respectively. Tilted boundary 
surfaces exhibit a larger entanglement per contacting site 
due to the increased connectivity. 

In the case of blocks separated by one site, we should 
use instead the full Eq. (|25|) with the second order ex- 
pression (j30j) . For parallel blocks with n sites at sep- 
aration s = 1, the n negative eigenvalues of the ma- 
trix (j2Tj) become, neglecting edge effects and setting 
a+ = |A+|/(4A), 

ff C K-(2aZo x -ol). (52) 

For blocks separated by one site through a 45° tilted 
surface of n modes, a discrete Fourier transform leads, 
neglecting edge effects, to (see Eq. (|C3|> ) 

/f ' C « -{2[a£„ + a 2 + a 2 + 2a x „(a ai + a,) cos 2s* 

+2a,a y cos^] 1 / 2 -cr 2 }, (53) 

where = a+a^, a x ,y = ^t^y + anu a k denotes 

the expression (|40p for m = n. In the isotropic case, Eq. 
([53]) becomes just 4<t(2o- + - a) cos 2 — . For the paral- 
lel and tilted subsystems c—d of Fig. tj] we then obtain, 
replacing sums by integrals and assuming a < 2a + , 

£0 ~ "cr(2cr+ - cr) , (54a) 
^ ~ 2ncr(2cr+ - cr) . (54b) 

Hence, the logarithmic negativity of the tilted case is, 
remarkably, twice that of parallel blocks when separated 
by one site, instead of 4/ir ~ 1.27 as in the contiguous 
case (Fig. |4]). Since they are a second order effect, Eqs. 
(f54| are not of the form (|45|) but rather 

£^ c w Lmcr(2cr + - a) , (55) 

if m is again the number of connections with the en- 
vironment per mode. They scale, therefore, with the 
measure \dA\2 (Fig. [5]). For larger separations s the neg- 
ativity vanishes at second order in A/A, as Fg C will be of 
higher order while the counter terms Gg and Gc remain 
of second order for sites at the surface. Consequently, 
the negativity becomes vanishingly small for s > 2. 




123456789 



FIG. 4. (Color online) Top: Exact and asymptotic logarith- 
mic negativities (Eqs. (|50l) - (|5ip ) for subsystems of the type of 
fig. [3] (for 10 x 10 blocks) as functions of A/A c . Tilted blocks 
exhibit a larger negativity per contacting site. Bottom: The 
tilted to parallel logarithmic negativity ratio for separations 
s — (a-b) and 1 (c-d). It is asymptotically 4/-7T in the 
contiguous case and 2 for one site separation. 

We finally remark that previous expressions are inde- 
pendent (at leading order) of the width d of the blocks 
(assumed finite), provided d > 2. In the case of two lines 
(d = 1, Fig. [S]), the extra interaction with the environ- 
ment at the other side of the line leads to an additional 
negative second order contribution in Eq. (|26[) . Hence, 
while it can be neglected in the case of contiguous lines, 
it will double the negative term in Eqs. ([54]) in the case 
of lines separated by one site, leading to the lower values 

~ na{2a + - 2cr) , (56a) 
I w 2na{2a+ - 2a) , (56b) 

which are now valid for a < a + . In this case negativity 
will vanish at second order if a > <j + . Edge effects in 
Eqs. ([54]) - ([56]) are also of second order and lead, using 
Eq. ([25]) . to a negative correction —2a 2 . 

All present expressions can be directly extended to 
three dimensions if present subsystems are extended 
parallel- wise along the z axis, replacing n by nn z . 

C. The fully connected case 

The evaluation of singular values is also straightfor- 
ward in the opposite case of a fully and uniformly con- 
nected system of n modes 0, [30l - [33 | (LMG type model 
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FIG. 5. (Color online) Top: Asymptotic and exact values 
of the logarithmic negativity of two contiguous n x n blocks 
with parallel and tilted boundary surfaces, scaled with |cL4|i 
in the main panel and \dA\2 in the inset, with OA = dBC\dC. 
Bottom: Same details for two blocks separated by one site. 
The appropriate scaling is verified to be \dA\i in the top panel 
and dA | 2 in the bottom panel. 
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FIG. 6. (Color online) Top: Parallel and tilted lines separated 
by one site. Due to the extra interaction with the environ- 
ment, the associated negativity (Eqs. ([56}) is lower than that 
of the corresponding blocks (c, d) of Fig. |3] (Eqs. |54j|), as 
appreciated in the bottom panels for the parallel case. 



where n A = n — n A , leading to 



where 



a; 



1.7 



(57) 



Here we can also compare with full exact results, since it 
is exactly and analytically solvable [H, . The present 
system can be used to describe entanglement between 
systems whose separation is small in comparison with 
the correlation length. 

In the present case, the matrices F^ are obviously con- 
stant for i ^ j, i.e., 



Ft 



(58) 



and the entanglement between disjoint subsystems B and 
C will just depend on the number of sites in B and in C, 
being independent of their separation or shape. 

The matrix Fg c will then have just a single non-zero 
singular value V disjoint B,C, namely (see Appendix) 

o BC = Vn&*\ F i\- (59) 



In the approximation (|19l) we then obtain a single non- 
zero symplectic eigenvalue for any global bipartition 



-\2 



£ 



A.A 



nAn^lFi rlog{n A n A (F 1 )7 



which corresponds to an area |<9*4|2 = n A n A in (|32[) (here 
Ma = 1 for i ^ j). 

Similarly, we obtain a single negative symplectic eigen- 
value for any pair of subsystems B,C, given by (|5"5|) or, 
in the complete approximation (f2T j) - (|25|) . by 



f' L w -Vn&*\ F i\ + l\Fx\\n B n e + n c n c ) , (61) 

with £g C = —fe,c- The second term in (|6T|) becomes 
important for small subsystems in a large environment 
(ng,nc <C n). Otherwise it can be neglected, in which 
case (|6TT) corresponds to \dA\i = \fn^h~c in (|33)) ~ (|36)) . 
For the scaling (|57p . F^ is proportional to so that 
Eqs. (f5T)| - ([6"Tj) remain finite for large n. The scaling is 
then again as in Eqs. pi)) - (P5|) with L = 1, j = and 
m = n A n A for global partitions or m = njsnc f° r a pair 
of subsystems. 

The previous picture is, remarkably, also that of the 
exact treatment, where there is a single non-zero sym- 
plectic eigenvalue f A for any subsystem A, given by 



f A = 



Fi n A n A /r 



1 

2 ' 



(62) 



(see ref. HH and Appendix [C]) . Here we have used the 
(60) local basis where Fq = in Eq. (|57p . in which case Eq. 
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© leads to 

Ff+F+/n={F-f. 

A first order expansion of (|62| in leads then to f A ps 
F\ n A n Al n i which coincides with Eq. (|60p since for weak 
coupling F^/n pa (F^) 2 . 

In the same way, the exact partial trans pos e of T>bc 
has a single negative symplectic eigenvalue [24j | 

(63) 

where /3 e ,c = n B n c /n, 7 BjC = + n c )(n — rig — 

n c)/ n + 2(3b,C (see appendix ICj) . Expansion of (|63|) up 
to first order in .F+ leads then exactly to Eq. (|61|) set- 
ting i 7 "^ ps n(i 7 ' 1 ~) 2 . The present approximate scheme 
allows then to immediately determine the weak coupling 
expressions (|60|) - (j6Tj) and to rapidly identify their be- 
havior with sizes n^, n B , nc- The exact value of F^ (to 
be inserted in (p)2"j) - (j6"3"j) ) is given in the appendix (Eq. 
(|U4p). Up to first order in A_ we obtain F-f ~ 2<n-i)x • 



IV. CONCLUSIONS 

We have shown how entanglement properties of weakly 
correlated gaussian states can be recast in terms of the 
singular values of a sub-block of the generalized con- 
traction matrix associated with the state. This allows 
to obtain in a quite simple way analytic expressions for 
both the entanglement entropy between complementary 
subsystems and the logarithmic negativity for non com- 
plementary subsystems, which imply distinct area laws 
for these two quantities in the case of short range or 
constant couplings. Several illustrative examples were 
considered, which show the dependence of these laws on 
the geometry, connectivity and separation between the 
subsystems. A final comment is that through applica- 
tion of the bosonic RPA formalism 0, H^, [33[ or other 
bosonization treatments, [H|,[3l|, the present scheme can 
be applied to weakly interacting spin systems. Moreover, 
it can in principle be also implemented in phases exhibit- 
ing symmetry-breaking at the mean field level (i.e., fields 
below the critical field in attractive XY or XYZ chains) 
away from the critical field, provided the proper multi- 
plicity corrections accounting for the different degenerate 
mean fields [24[ is taken into account. Such application 
is currently being investigated. 

The authors acknowledge support of CONICET 
(NC,JMM) and CIC (RR) of Argentina. 



Appendix A: Perturbative expansions for the 
symplectic eigenvalue problem 

In this work we have used perturbative results which 
are not necessarily trivial and which can be obtained fol- 



lowing techniques similar to those employed in the per- 
turbative diagonalization of the Dirac equation. We start 
with the symplectic diagonalization of the contraction 
matrix T>a of a subsystem A, which leads to the system 

FiU f -F^Vf = fU f , (Al) 
F^Uf-{l + F+)Vf = fV f (A2) 

where (^ ) is the symplectic eigenvector associated with 
the eigenvalue /. Eq. (IA2|) allows to write Vf as 

V f = [1(1 + .f) + F+}- 1 F A U f . (A3) 

Replacing (|A3|) in (|A1|) leads to the equivalent non-linear 
reduced diagonalization problem 

{F+ F A [1{1 + /) + F+]- l F A }U f = fU f . 

For small F , in agreement with the hypothesis that the 
state is weakly correlated, the symplectic eigenvalues / 
are small. Hence, at leading order Vf ~ F^Uf and we 
obtain the reduced standard eigenvalue equation 

(Fi-F A -F A )Uf = fU f , (A4) 

which leads to Eq. (18) and implies Eq. (21). If A is 
the whole system and the latter is assumed to be in the 
ground state of H, all / vanish and the relation F + w 
F~F~ (Eq. |H1)) is obtained. 

Let us now consider the Hamiltonian (|2|). The sym- 
plectic diagonalization of % entails the standard diago- 
nalization of A4H and leads to the system 

(A - A+)U U -A~V u =uU u (A5) 
&TU u -(A-A + )V u =uV u . (A6) 

In this case, HA - ^ is considered small. For a positive 
eigenvalue, the zero order approximation is obtained by 
neglecting all terms proportional to A - and V u , which 
are assumed small in comparison with U u , uj and A — A + . 
It leads to (A — A + )U U = wL^. which is a standard 
hermitian eigenvalue equation for U^. We then obtain 

V u = (A - A+ + uiyiR-Uu (A7) 
ps U*{Q + ujI)- 1 !! 1 A~ U u , (A8) 

(Eq. (J2SJ)), where we have written A - A+ ps UW\ 
with Q = diag(w a ) the diagonal matrix of eigenvalues. It 
should be noticed that if A is degenerate, A + will affect 
U considerably even if small. Ground state entanglement 
will remain however small since it depends on V. It can 
be also easily seen that expansion of Eq. (|A7|) up to sec- 
ond order in A* leads to Eq. (|30)) . 

Appendix B: Singular values 

The singular values a a of an arbitrary m x n matrix 
A are the square root of the non-zero eigenvalues of AA' 
or equivalcntly A 1 * A, which are both positive matrices 
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with the same non-zero eigenvalues. The singular value 
decomposition implies the existence of unitary matrices 
U, V such that A = UDV\ with D a diagonal ma- 
trix with diagonal elements a a or 0, and U, V unitary 
eigenvector matrices of AA* and A* A: AAXJ = UDD\ 
AAV = VD^D, i.e., AAU a = a 2 a U a , A AV a = a 2 a V a 
for the non-zero eigenvalues o~ a , with V a = A\J a ja a . 
For an hermitian A, a a = \X a \, with A Q the (non-zero) 
eigenvalues of A. 

The singular values determine the matrix m norm of 
A used in this work, defined as 



[Tr {A* A) 



m j 2 1 1 / m 



= C£< 



m\l/m 



(Bl) 



||j4||i is the trace norm, ||A||2 the standard Hilbert- 
Schmidt norm and ||^4||oo the spectral norm, which is 
just the largest singular value. 

The singular values a a of A also determine the non- 
zero eigenvalues of the hermitian (m+n) x (m+n) matrix 



B 



A 
A 



(B2) 



which are ±er Q , since B 2 = ( AA ^ A ) has eigenvalues a\. 
Eigenvalues ±c a correspond to normalized eigenvectors 
{ ± U vJ/V2, with AAU a = a 2 a U a , V a = AU a /o a and 
UlUp = 6 a p,V*Vp = 6 aP . 

These results first imply that the non-zero eigenval- 
ues of the matrix (fT5)l are the square of the singular 
values a A ' A of F~ A , as F A A = {F AA )\ implying 

a a ' A = °"a ■ They also entail that the negative eigen- 
values of the matrix (|21j) are minus the singular values 



a c,a f f when Gb and Gq are neglected. 



Appendix C: Evaluation of singular values 

In the first order approximation (|29[) . the matrix Fg C 
for first neighbor couplings and disjoint contiguous blocks 
B, C with n contacting sites, has elements of the form 



(FB. c )ij = fU-i) 



(CI) 



if the sites are adequately ordered, where f(l) = <5jo0^_l 
for parallel and f(l) = <j x 8iq + <J y 8n for tilted blocks. 



For blocks separated by one site, we should use the 
second order approximation (|30p , which leads again to a 
matrix of the form (|C1|) . with f(l) = 25ioo- + a for par- 
allel blocks and f(l) = 2[a^a x Sio + {°~i°~y + &y0x)8li + 
a+o-ySa] for tilted blocks. In all previous cases, F^ C F^ B 
is an hermitian matrix with elements of the form 



^2f(k-i)f(k-j)=g(i 



(C2) 



if edge effects are neglected, where g(l) = J2k f(k)f(k + 
I) = g(—l). Such matrix can then be exactly diagonalized 
(neglecting edge effects) by a discrete Fourier transform 
[32l |. leading to eigenvalues a 2 = J2 l g{l)e l27rkl / n , where 
k = 0, . . . , n— 1 and n is its dimension (this result is exact 
ifg(-l) = g(n—l)). For real g(l), as in the previous cases, 
we then obtain 



3(0) 



2 g(l) cos 
z>o 



2nk 



(C3) 



which leads to Eqs. ([iD |) -([g3 )l (in the case of C — A with 
A the tilted block, the final matrix F^jFj^j^ is again of 
the form JC2J). 

In the fully connected case, the exact singular values 
arise immediately as the matrix Fg C is just a rank 1 

c V i, j, which therefore has 



constant matrix, i.e., F, 



B,c 



Fr? n Fr> 



a unique non-zero singular value a = yjnjsnc \ c \ ■ ± b c ± C B 
is a ng x rank 1 matrix with constant elements ne|c| 2 , 
whose unique non-zero eigenvalue is ngnc \c\ 2 due to trace 
conservation. 

The full exact symplectic diagonalization can also 
be performed (see appendix in [24j for details). We 
quote here that the exact symplectic eigenvalues of 
the reduced state of L sites for the couplings ((57)) arc 

5 and (To = 



0"! 



1)2 
2> 



(F - 



LF,-) 2 



( F o + k) 2 - ( F o) 2 - \ ( L - 1 fold degenerate). For 
a pure global state, Co = 0. In the local basis where 
= 0, this implies F + = 0, which leads to Eq. (gg) . 



In the same way, we obtain Eq. . The exact value of 
the present F^~ was also evaluated in [24| in terms of a 
parameter A (Ff = A/(2n)): 



where ui 



F 



i(A 2 



2 ) 



4(n— 1)wqWi 



(C4) 



CJ 



= y/(\- A x )(\- A y ) and 



= V / ( A + ^)( A + ^r)' with A ± = ( A - ± A vV 2 - 
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